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ABSTRACT 


The  numerical  method  of  characteristics  is  applied  to  the  plane^ 
cylindricali  and  spherical  blast  wave  problems.  The  calculation  begins 
at  a  constar^t  time  line  within  the  blast  wave.  Along  this  line^  all 
three  flow  variables  are  prescribed  according  to  the  similarity 
solution.  Along  a  fixed  back  boundary,  which  lies  between  the  shock 
front  and  the  time-axis one  flow  variable  is  prescribed.  On  the 
other  boundary I  the  strong  shock  equations  apply. 

When  compared  to  the  exact  similarity  solution »  the  results  of 
the  method  of  characteristics  are  found  to  be  accurate  to  within  1%, 
for  all  variables  along  the  shock  fr.ont^  after  a  pressure  drop  of  99%. 
AlsO|  h^-type  extrapolation  of  calculated  results  always  improves  the 
accuracy f  whereas  the  h-type  extrapolation  may  not. 


V 


nomenclature 


c  •  sound  spaed 

E  ■  parameter  proportional  to  the  total  energy  within  the  wave 

h  ■  parameter  proportional  to  the  mesh  size 

k  ■  position  ratio  of  back  boundary 

P  ■  pressure 

t  ■  time  variable 

u  ■  particle  velocity 

U  »  shock  wave  velocity 

X  «  space  variable 

y  »  specific  heat  ratio 

V  ■  geometric  parameter  (v  *  1^  2,  or  3  corresponds  to  plane, 
cylindrical,  or  spherical  flow,  respectively) 
p  *  density 

•  constant  density  outside  wave  zone 

Subscripts 

s  *  at  the  shock  front 

b  «  at  the  back  boundary 


!.  INTRODUCTION 


In  the  hydrodynamic  regime  of  hypervelocity  impact  problems, 
the  governing  equations  are  the  conservation  equations  for  unsteady, 
compressible  fluid  flow.  Due  to  the  nonlinearity  of  these  equations, 
very  few  closed  form  analytical  solutions  exist  and,  for  most  problems, 
numerical  methods  are  the  only  resort.  Among  the  numerical  methods, 
the  two  most  commonly  used  are  the  finite-difference  method  and  the 
method  of  characteristics.  In  the  finite-difference  method,  an  arti¬ 
ficial  viscosity  must  be  added,  either  explicitely,  as  in  the  von  Netunan- 
Richtmyer  q  method^,  or  iroplicitely ,  as  in  the  Lax  method^.  By  the 
introduction  of  artificial  viscosity,  discontinuities  (shocks)  are 
spread  Into  narrow  regions,  across  which  flow  variables  change  rapidly 
but  continuously.  Flow  fields  with  a  complicated  shock  geometry  can 
be  calculated  without  extra  difficulty.  Because  of  the  smearing  of 
shocks,  however,  the  finite-difference  method  is  not  accurate,  especially 
in  regions  near  the  shock  front.  Also,  it  does  not  give  an  exact 
treatment  of  singularities  such  as  centered  rarefaction  waves.  As 
mentioned  in  Reference  3,  two  of  the  finite-difference  calculations  of 
hyperveJocity  impact  problems  gave  initial  peak  pressures  which  differed 
from  the  exact  values  by  at  least  15%,  In  the  method  of  characteristics, 
the  shock  is  calculated  by  the  exact  shock  equations;  therefore,  better 
accuracy  can  be  expected.  For  problems  with  simple  geometr)',  the 
method  of  characteristics  seems  to  be  more  suitable. 

For  solving  one-dimensional  unsteady  flow  problems,  the  method  of 
characteristics  as  described  in  textbooks  such  as  Courant  and  Friedrichs** , 
Shapiro^,  and  Rudinger^  has  been  used  widely.  Most  recent  accounts  of 
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this  method,  with  special  emphasis  on  high  speed  computation  by 
electronic  computers  may  be  found  in  books  by  Fox^  and  Hoskin®. 

Due  to  the  nonlinearity  of  the  equations,  an  exact  analysis  of  the 
error  involved  in  the  numerical  calculation  by  either  the  method  of 
characteristics  or  the  finite-difference  method  is  very  difficult. 

Makino  and  Shear®  showed  that  replacing  the  differential  characteristic 
equations  by  difference  equations  introduces  an  error  of  the  h^-type. 

But  the  type  of  error  in  the  variables  solved  from  these  difference 
characteristic  equations  is  still  not  known.  Without  an  exact  error 
analysis,  the  accuracy  of  any  numerical  method  can  be  estimated  by 
comparing  the  result  either  with  accurate  experimental  results  or  with 
an  exact  analytical  solution. 

In  the  present  paper,  calculations  of  the  blast  wave  problems  are 
made  using  the  numerical  method  of  characteristics,  and  the  results 
are  compared  with  the  exact  similarity  solutions.  If  the  region  in 
which  the  density  is  close  to  zero  is  excluded,  the  method  of  character- 
isitcs  yields  results  which  are  accurate  to  within  1%  on  the  shock 
front  when  a  99%  pressure  decrease  is  reached.  Different  types  of 
extrapolation  formulas  are  applied  to  three  sets  of  calculations  of 
the  same  problem,  each  set  having  a  different  mesh  size.  It  is  found 
that  the  h^-type  extrapolation  always  improves  the  results. 

These  results  of  the  blast  wave  calculations  indicate  that  for 
other  unsteady  flow  problems  which  do  not  include  extremely  low  density 
regions,  the  numerical  method  of  characteristics  can  be  expected  to  be 
of  the  same  accuracy. 

For  the  continuity  of  presentation,  the  basic  governing  equations 
and  the  corresponding  characteristic  equations  are  given  in  section  II; 
different  similarity  solutions  for  blast  waves  are  discussed  in 
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section  II1<  Discussions  on  the  forimilation  of  the  problem^  numerical 
calculations^  and  method  of  extrapolation  are  then  presented  in  the 
second  half  of  the  paper.  Appendices  A,  and  C  are  included  mainly 
for  easy  references  although  some  of  the  material  represents  new  ways 
of  presenting  these  topics.  Appendix  D  includes  an  8-digit  table  of 
properties  within  the  blast  waves «  This  table  may  be  used  to  evaluate 
the  accuracy  of  other  numerical  methods. 


; 

I I .  GOVERNING  EQUATIONS 


When  the  unsteady  flow  of  an  ideal  gas  depends  on  only  one  space 
coordinate,  the  governing  equations  in  Eulerian  coordinates  for 
continuous  flow  without  friction  and  heat  transfer  are  the  continuity 
equation 


(1) 


the  momentum  equation 
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and  the  energy  equation 
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where  v  has  a  value  of  1,  2,  or  3  corresponding  to  plane,  cylindrical, 
or  spherical  flows,  respectively.  From  this  system  of  hyperbolic 
equations,  the  following  three  physical  characteristic  equations  may 
be  obtained  (see  Appendix  A). 

Along  I  and  II  dx 
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where  I  and  II  are  the  right-traveling  and  left  traveling  waves, 
respectively,  and  III  is  the  particle  path  line*  The  corresponding 
state  characteristics  are 

along  I  and  II,  ^  ^  ^ 

along  in,  dp»i(dp)  (7) 

In  these  characteristic  equations,  the  following  expression  for  sound 
speed  as  a  function  of  pressure  and  density  has  been  introduced: 

(8) 

After  eliminating  c^  by  using  eq„  (8),  eq»  (7)  may  be  integrated,  and 
the  result  indicates  that  p/p^  ,  or  the  entropy,  remains  constant  along 
each  particle  path*  However,  this  constant  is  in  general  different 
for  each  particle  path.  Notice  that  this  assumption  is  incorporated 
in  the  energy  equation  (3),  It  is  convenient  to  eliminate  the  density  p 
from  the  characteristic  equations  by  using  equation  (8),  The  state 
characteristics  are  then 

along  I  and  II,  du  •  ♦  —  (dp)  *  (v-l)  —  (dt)  (9) 

Yp  X 

2Y 

along  III,  12- *  (—)y-T  #  (10) 

^o 

The  constants  c^  and  p^  must  be  evaluated  for  each  particle  path. 
Equations  (4),  (5),  (9),  and  (10)  are  the  basic  equations  which  are  to 
be  used  in  the  method  of  characteristicSw 

Only  shock  waves  of  very  strong  intensity  will  be  used  in  our 
calculation.  'Fhis  restriction  is  consistant  with  the  similarity  solution 
of  blast  waves,  the  results  of  which  will  be  compared  with  our  method  of 
characteristics  solution.  The  Rankine-Hugoniot  relations  for  a  strong 
shock  moving  into  a  uniform,  stagnant  region  are 
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In  the  numerical  procedure  of  the  method  of  characteristics^ 
all  governing  equations  are  put  into  finite-difference  form.  Thus, 
equations  (4),  (5),  (9),  and  (10)  become 


along  I  and  II, 


along  III, 
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where  A  represents  the  (central)  difference  between  two  adjacent  points 
along  a  characteristic  and  the  barred  values  represent  the  average 
between  the  same  two  points. 


III,  SIMILARin  SOLUTIONS  OF  BLAST  WAVES 

The  self-similar  solutions  of  blast  waves  have  been  analyzed  by 
many  investigators  and  are  among  the  few  known  exact  solutions  of  the 
one-dimensional,  non-isentropic  (entropy  changes  from  particle  to 
particle),  unsteady  flows.  This  type  of  flow  will  be  calculated  by 
the  numerical  method  of  characteristics  and  the  results  compared  with 
the  analytical  solutions. 
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Taylor^ ^  solved  the  blast  wave  problem  which  consists  of  the 
instantaneous  release  of  a  finite  anotmt  of  energy  at  a  point  within 
a  stagnant^  ideal  gas  field.  The  subsequent  flow  was  non-isentropic^ 
with  spherical  symmetry  (u  •  0  at  x  ■  0),  and  with  a  strong  shock 
propagating  into  the  stagnant  gas.  By  utilizing  the  self«similar 
property  of  the  flow,  the  governing  partial  differential  equations 
were  reduced  to  ordinary  differential  equations  with  one  independent 
variable.  These  equations  were  then  integrated  numerically  and  the 
results  presented  in  tabular  form.  Lin^^  followed  essentially  the 
same  technique  and  solved  the  problem  of  the  instantaneous  release  of 
energ)*  along  a  line  and  the  ensuing  cylindrical ly  symmetric  flow. 

His  results  in  tabular  form  were  also  obtained  by  numerical  integration. 
Recently^  Lee^^  presented  a  unified  approach  which  encompasses  a  large 
number  of  similarity  solutions.  He  did  not  attempt  to  integrate  the 
resulting  ordinary  differential  equations. 

Sedov^^  and  Rouse^**,  among  others^  have  successfully  integrated 
the  ordinary  differential  equations  of  the  blast  wave  problems  and 
obtained  closed  form  solutions »  They  also  included  in  their  solutions 
the  case  of  plane  motion ^  in  addition  to  the  cylindrical  and  spherical 
cases.  Since  these  solutions  are  in  closed  form^  their  numerical  value 
can  be  calculated  to  any  desired  degree  of  accuracy  without  too  much 
difficulty.  On  the  other  hand,  the  solutions  of  Taylor  and  Lin  involve 
numerical  integration;  therefore,  a  high  degree  of  accuracy  is  difficult 
to  achieve.  Figure  1  gives  a  comparison  of  the  results  of  our  cal¬ 
culation  of  Sedov's  spherical  wave  solution  and  the  solution  obtained 
by  Taylor.  It  can  be  seen  that,  although  Taylor's  values  may  be 
satisfactory  for  other  purposes,  they  are  not  accurate  enough  to  serve 
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as  a  standard  to  evaluate  the  accuracy  of  other  numerical  Methods. 

Figure  10  of  Reference  14  gives  a  sinilar  cooparison  of  the  analytical 
solution  of  Reuse  for  the  cylindrical  case  and  the  nuaerically  integrated 
solution  of  Lin. 

It  is  evident,  then,  that  only  accurately  evaluated  closed  fore 
solutions  can  be  used  as  a  standard  to  gauge  the  accuracy  of  nuaerical 
Methods.  Gerber  and  Bartos^^  presented  a  nuaerical  table  calculated 
from  Rouse's  equations  for  cylindrical  waves  w5th  y  ■  7/S  and  5/3. 
However,  their  table  contains  too  few  data  points  for  our  purposes. 

For  easy  reference,  we  have  made  numerical  calculations  of  Sedov's 
solution  for  the  plane,  cylindrical,  and  spherical  cases  with  two  values 
of  Yi  7/5  and  5/3.  These  numerical  results  are  tabulated  in  Appendix  D, 
table  II.  Eight  significant  figures  were  used  in  all  calculations.  A 
comparison  with  calculations  made  for  a  few  points  using  twenty-two 
digits  indicated  that  the  data  in  Appendix  D  are  accurate  at  least  to 
the  sixth  digit.  Sedov's  solution,  for  y  »  7/5,  is  as  follows, 
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Bid  V  is  a  paraaeter  with  the  range  of  variation 
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Hie  properties  at  the  shock  front ,  which  are  indicated  by  the  subscript 


s  in  the  above  equations ^  are  given  by 
“s  *  ttW  ^57^^  *2"^ 
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The  position  of  the  shock  front,  at  any  time  t, 
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IV.  FORMJUTION  OF  THE  PROBLEM  FOR  NUMERICAL  CALCULATION 

In  iq)plying  the  numerical  method  of  characteristics  to  the  blast 
wave  problems,  it  is  desirable  to  avoid  the  origin  (x  «  0,  t  «  0)  and 
the  time  axis  (x  ■  0).  The  origin  is  a  singular  point  with  unbounded 
values  of  pressure,  particle  velocity,  and  sound  speed  and,  along  the 
time  axis,  the  sound  speed  is  infinite.  Only  a  region  of  the  flow 
field  which  does  not  include  the  time  axis  will  be  calculated.  This 
region  is  b.">unded  by  a  constant  (non-zero)  time  line  on  the  bottom, 
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the  shock  front  on  the  rights  and  a  back  boundary  curve  on  the  left 
trith  the  equation  ^  ^  \  »  where  k  is  a  constant.  (See  Figure  2.) 

To  pose  a  problea  which  is  Mtheutically  correct,  the  following 
boundary  and  initial  conditions  mat  be  prescribed  on  the  back 
boundary,  shock  front,  and  constant  tiae  line. 

The  constant  time  line,  t  ■  t^  in  Fig.  2,  is  considered  as  the 
initial  value  curve;  therefore,  values  of  all  three  dependent  variables 
p,  u,  and  c  Bust  be  specified  on  it.  These  values  were  calculated 
froB  the  exact  solution,  equations  (18)  to  (25)  and  equation  (8),  and 
considered  as  the  specified  values  for  the  solution  by  the  aethod  of 
characteristics.  The  back  boundary  curve  is  a  "tiae-like"  arc, 
according  to  the  terminology  used  by  Courant  and  Friedrichs**,  since 
two  of  the  three  characteristics,  II  and  III,  reach  this  curve  froB 
inside  the  region  R.  The  correct  boundary  condition  requires  that 
one  of  the  three  dependent  variables  be  specified  on  this  curve.  In 
our  case,  the  particle  velocity  u  was  calculated  froB  the  exact  solution 
and  used  as  the  one  specified  dependent  variable. 

The  shock  front  is  also  a  time-like  arc,  but  in  this  case  only 
the  I-characteristic  reaches  a  point  on  it  from  inside  region  R.  Also, 
there  is  no  prior  information  about  its  exact  position.  There  are  four 
unknown  quantities  to  be  determined  at  the  shock  front,  U,  p,  c,  and 
u.  Therefore,  the  strong  shock  conditions,  equations  (11),  (12),  and 
(13),  are  the  correct  boundary  conditions  at  the  shock  front  since 
these  three  equations,  together  with  the  equation  of  the  I-characteristic 
serve  to  determine  the  four  unknown  quantities.  The  position  of  the 
shock  path  is  determined  by  integrating  the  shock  velocity  U. 
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Ve  NUMERICAL  CALCULATION 


In  the  nunerical  calculation i,  a  nuober  of  equally  spaced  points 
on  the  constant  tiae  linc^  t  «  arc  chosen ..  and  the  properties 
at  these  points  are  calculated  fron  the  exact  solution o  A  network 
is  then  constructed  in  the  physical  plane  consisting  of  the  1  and  II 
characteristics  passing  through  these  points  and  the  reflections  of 
these  characteristics  from  the  back  boundary  and  shock  front 
Properties  at  the  mesh  points  of  this  network  will  be  determined. 

Since  the  equations  for  the  state  characteristics  involve  the  time 
variable  t,  and  the  physical  characteristics  involve  the  state 
variables  u  and  Cp  neither  the  physical  nor  the  state  plane  can  be 
independently  constructedc  An  iteration  scheme  must  be  adopted  to 
determine  a  mesh  point  in  both  the  physical  and  state  planeSr  i.e  ^  for 
each  point  the  five  unknowns  x  t..  p.  u.  and  c  must  be  solved  for 
simultaneously  by  some  iteration  process  The  iteration  scheme  adopted 
here  is  described  in  detail  in  Appendix  It  follows  essentially  the 
standard  procedure  given  in  textbooks  such  as  References  5  and  16. 

In  brief,  it  involves  the  construction  of  the  1  and  11  characteristics 
from  two  neighboring  points „  locating  their  intersection,  and  tracing 
back  a  path  line  until  it  interscvts  the  straight  line  connecting  the 
original  two  points «  At  this  intersection  point  the  entropy  is  calcu¬ 
lated  by  a  linear  interpolation  of  the  values  at  the  two  known  points 
The  calculated  results  for  plane ^  cylindrical,  and  spherical  waves 
are  presented  in  Table  I,  The  same  values  of  and  E  were  used  for 
all  three  cases.  For  the  factor  k,  which  determines  the  position  of 
the  back  boundary,  three  different  values,  0^5,  0,7,  and  0,8  were  used 
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For  all  wave  geoaetries  pi  «  ,0024497  slugs/ft  ,  E  *  1.5482649  x  10  ft-lb/ft 


for  the  plane^  cylindrical »  and  spherical  cases ^  respectively.  As 

the  value  of  k  decreasesi  i.e«|  as  a  point  moves  towards  the  time  axis^ 

the  value  of  the  density  p  approaches  zero  very  rapidly  and  the  sound 

speed  c  increases  rapidly  without  bounds  as  shown  in  Fig*  3*  Because 

of  these  extreme  values «  the  accuracy  of  the  numerical  solution 

becomes  inadequate  as  k  becomes  smaller*  The  mesh  of  the  characteristic 

network  in  the  physical  plane  is  of  irregular  size  when  these  extreme 

values  of  c  are  included.  In  Figure  4  a  coarse  network  for  the  case 

of  V  ■  k  >  0*6  is  shown*  For  practical  protlems  such  as  those 

encountered  in  hypervelocity  impacts,  exploding  wires,  or  atmoshpheric 

explosions,  the  initial  special  distribution  of  c  at  any  two  points 

seldom  differs  by  a  ratio  of  more  than  2.5,  and  the  initial  density 

is  seldom  close  to  zero.  Therefore,  the  value  of  k  was  so  chosen  that 

the  ratio  of  the  sound  speed  at  the  back  boundary  over  that  at  the  shock 

front  was  approximately  2.5  for  all  three  cases.  This  results  in 

the  aforementioned  values  of  0*5,  0*7,  and  0.8  for  the  plane,  cylindrical, 

and  spherical  cases,  respectively,  as  indicated  in  Fig.  3.  The 

initial  position  of  the  shock  front  was  aribtrarily  taken  as  2  feet 

for  all  three  cases:  therefore  a  different  value  of  t  was  obtained  for 

o 

each  case* 

The  numerical  procedure  was  programmed  in  Fortran  II  language  and 
the  calculation  performed  on  an  IBM  7040  computer.  The  maximum  com¬ 
putation  time  was  approximately  two  hours,  which  gave  a  pressure  at 
the  last  point  on  the  shock  front  of  about  1%  of  the  initial  pressure 
on  the  shock. 

The  accuracy  of  these  calculations  is  shown  in  Figures  5,  6,  and 
7*  Figure  5(a)  shows  the  physical  plane  of  the  plane  wave,  indicating 
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the  position  of  the  shock  fronts  the  back  boundary ^  and  the  last  line 
of  nuaerical  calculation.  The  percent  relative  error  of  p»  u,  and  c 
at  points  on  the  last  line  of  calculation  as  compared  to  the  exact 
solution  is  plotted  in  Fig.  5(b).  The  percent  error  of  t»  p»  u^  and  c 
along  the  shock  front  is  shoirn  in  Fig«  S(c).  These  errors  along  the 
shock  front  were  obtained  by  coaparing  the  calculated  values  of  these 
variables  with  those  of  the  exact  solution  at  the  saae  x  location. 
Since  the  location  of  the  numerically  calculated  shock  front  deviates 
slightly  from  the  exact  shock  fronts,  the  properties  are  not  co^)ared 
at  the  same  value  of  time;  however ^  this  error  in  time  is  quite  small. 
The  maximum  error  along  the  shock  front  for  all  variables  is  less  than 
.6%;  while  on  the  last  line  of  calculation  the  maximum  error  is  less 
that  3%.  Similar  results  for  the  cylindrical  and  spherical  cases 
are  shown  in  Figures  6  and  7^  respectively. 


VI.  EXTRAPOUTION  OF  NUMERICAL  RESULTS 

It  is  shown  in  Appendix  C  that  the  error  introduced  by  replacing 
the  differential  equations  (4) «  (5)^  and  (9)  by  finite-difference 
equations  (14),  (IS),  and  (16)  is  of  the  h^^type,  where  h  is  a  para¬ 
meter  proportional  to  the  mesh  size.  Also,  the  linear  interpolation 
process,  which  is  used  in  the  solution  of  the  difference  equations, 
is  known  to  be  correct  to  the  order  of  h^.  The  boundary  conditions 
and  initial  conditions  used  in  our  problem  are  either  prescribed 
functions  or  exact  strong  shock  conditions.  Since  all  errors  introduced 
by  the  numerical  solution  are  of  the  order  h^,  it  seems  reasonable  to 
assume  that  each  of  the  numerically  calculated  variables  has  an  error 
which  is  of  the  h^-type,  i.e., 

-  p  ■  ♦  ♦jh®  ♦  ...  (26) 
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where  p  is  the  exact  value  of  pressure,  p  is  the  approximate  pressure 
at  the  same  point  calculated  with  a  mesh  size  h,  and  are  quantities 
which  are  independent  of  the  mesh  size  and  depend  only  upon  the  location 
of  the  point  in  question.  Equations  similar  to  equation  (26)  are  also 
assumed  to  hold  for  c  and  u. 

With  this  order  of  error  assumed,  it  is  possible  to  apply 
Richardson’s  extrapolation  method to  the  calculated  values  at  any 
point  with  two  different  mesh  sizes  and  obtain  new  values  which  are 
more  accurate.  If  two  values  of  the  pressure  at  a  point  are  calculated, 
p^  corresponding  to  a  mesh  size  h^,  and  p^  corresponding  to  a  mesh  size 
h^,  then  the  extrapolated  value  may  be  determined  from  the  following 
equations  which  are  obtained  by  truncatiiig  equation  (26) , 


p.  •  Pi  *  *i‘'i 

P.  ’  P2  *  *1^1 

Elimination  of  from  these  equations  gives 


(27) 


-h2 


(28) 


In  this  equation,  p  represents  the  extrapolated  value  of  p.  Similarly, 

c 

if  calculations  with  three  different  mesh  sizes,  h^  h2,  and  h3,  arc 

performed,  the  following  equations  may  be  used  to  solve  for  p  , 

e 


Pe  -  P,  ♦ 
Pe  ”  Pj  * 

P,  -  P3  ♦ 
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It  must  be  mentioned  thatp  although  the  finite-difference  equations 
Involve  errors  of  the  h^-type,  there  is  no  proof  to  show  that  the  errors 
in  the  calculated  values  of  p^  u,  and  c  themself  are  also  of  the  h^-type. 
Equation  (26)  is  merely  an  assumption.  However^  for  the  present 
problem^  as  demonstrated  in  the  numerical  examples  below^  the  h^-type 
extrapolation  does  improve  the  numerical  results. 

Roberts  used  an  extrapolation  formula  which  is  based  on  the 
assumption  that  all  errors  in  numerical  solutions  are  of  the  h-type; 
or,  in  our  notation ,  his  assumption  becomes 

where  are  functions  independent  of  h  and  dependent  only  upon  the 
position  of  the  point  in  question.  The  resulting  extrapolation  formula 
for  two  sets  of  calculations  is 


Pe  .  hj 


(31) 


For  co«q)arison  purposes^  both  h^-type  and  h-type  extrapolations 
are  applied  to  the  present  numerical  results  for  the  plane  wave  case. 

In  addition  to  the  101  initial  point  (100  initial  segments)  calculation » 
two  more  sets  of  calculations,  using  50  and  25  initial  segments, 
respectively,  were  performed.  The  mesh  si;^es  of  these  three  cases 
are  designated  as  h  ■  1/4,  h  «  1/2,  and  h  •  1,  respectively.  The 
errors  in  p  along  the  shock  front  for  these  three  cases  are  shown  in 
Fig.  8.  It  can  be  seen  that  the  error  decreases  as  the  mesh  size 
decreases.  Several  sets  of  extrapolated  values  of  the  pressure  along 
the  shock  front  were  obtained  by  using  different  combinations  of  the 
three  sets  of  calculations  with  different  mesh  sizes.  These  extrapolated 
values  of  p  were  compared  with  the  exact  values  and  the  error  plotted  in 
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Figure  8.  The  curve  labeled  h^(l,  1/2)  represents  the  error  in  the 
extrapolated  valai  of  p  obtained  froai  the  h^-type  formula,  eq«  (28), 
by  using  the  pressures  from  the  corresponding  calculations  with  mesh 
sizes  h  ■  1  and  h  ■  1/2  ;  h(l,  1/2,  1/4)  represents  h-type  extrapolation 
of  pressures  corresponding  to  mesh  sizes  h  ■  1,  1/2,  and  1/4,  etc*  The 
h(l/2,  1/4) extrav^lation  gives  the  least  error,  while  the  other  three 
h«extrapolations  give  errors  equal  to  or  larger  than  the  original 
calculation  of  mesh  size  h  •  l/4«  On  the  other  hand,  all  four  h^* 
extrapolations  give  errors  which  are  smaller  than  the  h  «  1/4  case. 

In  other  words,  for  the  present  problem,  the  h^^extrapolations  always 
improve  the  results,  whereas  the  h-extrapolations  do  not  necessarily 
iflq3rove  the  results.  The  results  also  indicate  that  extrapolations 
of  three  sets  of  calculations  are  not  necessarily  better  than 
extrapolations  of  two  sets  of  calculations.  Extrapolation  curves 
for  the  variables  u  and  c  give  essentially  the  sane  results  as  those 
of  pressure  and,  therefore,  are  not  presented. 

VII.  CONCLUSION 

The  numerical  method  of  characteristics  produces  very  accurate 
results  when  applied  to  blast  wave  problems  if  the  region  where  the 
density  is  close  to  zero  is  excluded.  For  other  one-dimensional 
unsteady  fluid  flow  problems  which  do  not  include  zero  density  regions, 
it  is  reasonable  to  expect  that  the  numerical  method  of  characteristics 
will  also  produce  accurate  results. 

According  to  our  numerical  calculations,  h^-type  extrapolation 
always  improves  the  accuracy,  whereas  h-type  extrapolation  may  not. 
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For  the  purpose  of  evaluating  the  accuracy  of  any  nuaerical 
method  used  for  solving  unsteady ^  non*isentropic  flo^s,  the  similarity 
solution  of  blast  waves  may  be  used  as  a  standard  for  comparison. 

Since  the  similarity  solutions  are  quite  lengthy^  they  have  been 
computed  for  two  values  of  and  the  numerical  results  are  presented 
in  Appendix  D  for  easy  reference. 
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Figure  1  Comparison  of  Taylor’s  solution  (numerical) 
with  Sedov’s  soltuion  (analytical)  for  the 


spherical  blast  wave,  -  -  -  .  -  Taylor 
20  ■"  Sedov 


I 

I 


Figure  2  Physical  Plane  Showing  Region  of  Numerical  Solution  for  Blast  X 

Nave  Problems;  on  the  back  boundary  (x^  ■  kx^) »  u  is  prescribed; 
on  the  constant  time  line»  p»  u,  and  c  are  prescribed;  on  the 
shock  fronti  the  strong  shock  equations  (11)  to  (13)  apply. 
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RELATIVE  ERROR  ,  %  TIME  I ,  t«c 


Figure  5- Plane  Wove  Error  Anolytie 

....  ,  ¥  ‘ 


RELATIVE  ERROR  ,  %  TIME  t  .  ttc  x  10 


Figure  6-*  Cylindrical  Wave  Error  Analysis 
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RELATIVE  ERROR.  %  TIME  4  . 


Figure  7—  Spherical  Wave  Error  Analysis 
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RELATIVE  ERROR 


h(l  i) 


Figure  8  Curves  of  Relative  Error  in  Pressure  for  tho  Plane  Blast  Wave 
(k  •  .5) 


- - results  of  numerical  calculation  using  grid 

sizes  h  -  1,  h  -  1/2,  and  h  -  1/4 
—  —  —  —  —  —  —  h-type  extrapolation 
■  '■  h^“typc  extrapolation 

h(l|  1/4)  indicates  h-type  extrapolation  of  the  two  results  obtained 
from  grid  sizes  h  «  1.  and  h  ■  1/4. 
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APPENDIX  A 


Derivation  of  Characteristic  Equations 

In  this  appendix 9  the  application  of  the  method  of  characteristics 
to  the  one*dinensional  flow  of  an  ideal  gas  with  area  change,  friction, 
and  heat  transfer^  is  presented.  The  governing  equations  are  the 
continuity  equation 

the  momentum  equation 

and  the  energy  equation 

where  F  is  the  wall-friction  tern  and  q  represents  the  heat  transfer 
rate*  From  the  ideal  gas  equation  of  state  and  tne  definition  of 
sound  velocity,  we  have 

C  -  (A, 4) 

With  this  equation,  the  properties  a,  p,  and  p  may  be  treated  as  dependent 
variables  which  are  governed  by  equations  (A.l),  (A*2),  and  (A, 3).  For 
regions  where  these  variables  are  continuous,  three  equations  for  the 
total  differentials  du,  dp,  and  dp  may  be  written.  Combining  these 
differential  relations  with  the  governing  equations,  (A.l)  to  (A. 3),  we 
obtain  the  following  set  of  six  equations  which  may  be  solved  for  the 
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six  partial  derivatives  of  p,  and  p« 


+dx|U+<jtf 

+  <ixf+(ltf 

+<lt^ 

f  ^  +  if 

1+  f f 


dU 

dP 

^  ,..s, 

-F 

(v-miiF+<f 


Solving  these  equations  for  au/dx  by  Cramer's  RulCi  we  obtain 

»=# 


where 

dx  dt  0  0  0  0 

0  0  dx  dt  0  0 

0  0  0  0  dx  dt 

pOOOul 
u  1  1/p  0  0  0 

0  0  u  1  -uc^  -c^ 

and 


‘(liH-MiLOdt-d  jJBii-Cldt  -dx] 


du  dt  0  0  0  0 

dp  0  dx  dt  0  ^ 

dp  0  0  0  dx  dt 

•pu/A  dA/dx  0  0  0  u  1 

-F  1  1/p  0  0  0 

(Y-l)p(uF^q)  0  u  1  -uc^  -c^ 


=(a<lt  +<(U)fU-dt  -dx) 
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The  vanishing  of  D  produces  the  three  physical  characteristics 


I-characteristics 

^ru-cc 

(A.-) 

Il-diaracteristics 

^=tL-C 

(A.  8) 

Ill-characteristics 

II 

{= 

(A.9) 

Notice  that  the  Ill-characteristics  are  also  path  lines.  Along 
these  physical  characteristics ^  the  derivative  du/dx  is  indeterminate 
and  therefore  may  be  discontinuous.  To  insure  that  du/dx  is  indeter- 
and  not  infinitOi  the  numerator^  N,  of  equation  (A.6)  must  also  vanish 
along  the  characteristics.  Along  the  Ill-characteristic^  (udt-dx)  *  0^ 
and  N  vanishes  identically.  Along  the  I  and  I I- characteristic » 
dx/dt  •  u  i  C|  and  the  vanishing  of  N  produces  the  following  state 
characteristics  (or  coiiq)atibility  equations). 

(a.io) 

The  upper  signs  refer  to  I-characteristic$  and  the  lower  signs  refer 
to  ll-characteristics. 

When  equations  (A, 5)  are  solved  for  du/dt,  dp/dx,  and  dp/dt  , 
the  vanishing  of  the  numerators  yields  results  which  are  identical  to 
those  for  du/dx  .  When  equations  (A.S)  are  solved  for  dp/dx  and  dp/dt  , 
the  numerators  of  these  solutions  do  not  contain  the  common  factor 
(udt-dx).  The  vanishing  of  these  numerators  yields ^  in  addition  to 
equation  (A.IO),  the  Ill-state  characteristic, 

d /=-^[dP-(y-Jl/(UF+(jUt]  (A. 11) 
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The  solutions  of  equations  (A,5)  for  the  other  five  derivatives 
are  given  belov  for  reference. 

^  -^U  dt-<lx){[(y-lVuf  dxdt  +[F(u.dt-dx) 

-  lidu  -^dx+(u'-c’)  du  dt} 
■^=>()t-dx)lInr-ivruF*(i)dt-^’^dt-dij(udt-<W 

+/c'dtrFdt+du)} 

^  ^{udt-dxijfSj^'^dx-ry-iWurt^idx+udpjx  '*•”> 

(udt-dx)-A^Fdxdt-c*dPdt->’c*diidx]' 
■^=-^(adt-dx){(--^f‘^dt-d/)(udt-dx)-'/F((WVdudt} 
+(X-ll>’(ttF +i()(dtf+  d/(()tPc*  -dPCdti*) 

dx +U<l/^(tldt  -dx)  ->>dx(  F+d  u)} 

-  (i’-ll/’fil  F+</)  dx(dt)=^-d  Pdxdt  -  iic’d/Wt)’) 


For  frictionless  flow  without  heat  transfer »  both  F  and  q  vanish. 


The  quantity  (l/A)(dA/dx)  can  be 

evaluated  for  the 

following  specific 

types  of 

one-dimensional  motion: 

1) 

Plane  motion 

A  *  coastant, 

1  dA  « 

A33r“  ^ 

2) 

Cylindrical  motion 

A  «  2Tix-t  , 

1  dA  1 

X33r‘7 

3) 

Spherical  motion 

A  ■  4ttx^  , 

i  .  i 

A  3x  “  7 

1 

I 

a 

I 

I 
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where  x  represents  the  redial  coordinate  for  the  cylindrical  and 
spherical  cases.  These  results  nay  be  sumnarized  by  writing 

where  v  «  1^  2,  or  3  corresponds  to  plane ^  cylindrical ^  or  spherical 
iBOtionSf  respectively.  The  physical  and  state  characteristic 
equations  for  one-dinensional  flow  of  an  ideal  gas  without  friction 
and  heat  transfer  become 

Physical  Qiaracteristics  State  Characteristics 

{(Wjj  W  P)j|[ 

It  is  convenient  to  eliminate  the  density  p  from  the  characteristic 
equations  by  using  equation  (A.4).  When  this  substitution  is  made^ 
and  the  third  state  characteristic  equation  integrated^  the  following 
characteristic  equations ^  in  terms  of  and  c^  are  obtained. 

Physical  Characteristics  State  Characteristics 


The  constants  c  and  p  must  be  evaluated  for  each  Ill-characteristic. 
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APPENDIX  B 


Iteration  Procedure 

In  this  appendix^  the  solU(.ion  of  the  finite-difference 
characteristic  equations  (14)  to  (17) ^  for  the  boundary  conditions 
specified  in  Section  IV,  is  described. 

The  basic  step  in  the  solution  of  unsteady  flow  problems  by  the 
method  of  characteristics  is  the  determination  of  all  variables,  u, 
c,  p,  X,  and  t,  at  a  new  point  when  these  variables  are  known  at  two 
existing  points.  Repi  ned  application  of  this  step,  using  suitable 
initial  and  boundar  conditior^s,  produces  the  solution  of  the  problem. 
For  the  blast  wave  problem,  three  different  types  of  points  exist, 
each  requiring  a  different  procedure  in  the  basic  step.  These  three 
types  of  points  are 

1)  interior  points,  which  lie  within  the  region  of  numerical 
solution, 

2)  points  on  the  shock  front,  and 

3)  points  on  the  back  boundary. 

The  three  procedures,  which  were  used  to  solve  for  the  three  types  of 
points,  are  outlined  below. 

1)  Interior  Points 
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Referring  to  Figure  all  variables  are  assumed  knoim  at 
points  1  and  2,  and  the  variables  at  point  3  are  to  be  determined. 
Point  A  is  located ^  in  the  physical  plane ,  at  the  intersection  of 
the  111-characteristic  (path  line)  passing  throu^  point  3  and  the 
straii^t  line  joining  points  1  and  2.  The  properties  c^^  and  p^, 
at  point  A,  are  determined  by  a  linear  interpolation  of  the  properties 
at  points  1  and  2.  The  finite-difference  forms  of  the  physical 
characteristic  eqiiations  are,  from  equations  (14)  to  (17) 


along  1 

(B.l) 

along  II 

iio*  c„ 

(B.2) 

along  III 

it 

II 

(B.3) 

The  corresponding  state  characteristics  are 

along  I 

Uj-Ur  - 

(B.4) 

along  II 

U3-Ua= 

(B.S) 

along  III 

(B.6) 

where  the  barred  quantities  represent  average  values,  i.e., 


The  equation  of  the  line  between  points  1  and  2  is 


and  the  linear  interpolation  formulas «  which  are  used  to  determine  the 
properties  at  point  are 


(B.8) 

(B.9) 

(B.IO) 

From  these  ten  equations^  the  ten  unknowns^  x  .  t  ,  u  ^  c  ,  p  «  x., 

3  3  3  3  3  A 

^A*  ^A*  ^A*  determined.  The  solution  is  accomplished 

by  the  following  iteration  procedure. 

First  Iteration 

A  single  prime  is  used  to  designate  the  first  approximation  of 
the  variables. 

Equations  (B.4)  and  (B.S)  are  solved  simultaneously  for  the 
properties  u^  and  p^  by  neglecting  terms  containing  t  and  using  the 
properties  at  point  1  to  replace  the  average  properties  between  1  and 
3,  and  the  properties  at  2  to  replace  the  average  properties  between 
2  and  3.  These  values  are  then  designated  by  u  '  and  p  An  estimated 
value  of  c^,  which  is  called  c^^  is  calculated  from  the  equation 


Nextf  equations  (B.l)  and  (B.2)  are  solved  simultaneously  for  x  '  and 
t^'  by  using  the  previously  determined  values  of  c^  and  u^'.  Equations 
(B.3)  and  (B.7)  are  then  solved  simultaneously  for  x^'  and  t^'  by 
using  u^’  in  place  of  u^^  .  The  properties  u^*,  c^’ ,  and  p^*  are 
determined  from  the  interpolation  formulas,  equations  (B.8),  (B.9),  and 
(B.IO),  and  the  calculated  vaijw  of  x^'.  The  quantity  is  calculated 
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by  inserting  the  values  and  p^*  into  eq  lation  (B.6).  The 

first  approxifluition  of  the  properties  at  points  3  and  A  is  now  complete. 


Second  Iteration 

Once  the  first  approximation  is  established,  the  following 
procedure  may  be  used  to  obtain  the  second  and  subsequent  approximations 
of  the  variables  at  point  3.  A  double  prime  is  used  to  designate  the 
second  approximation. 

The  values  of  u  **  and  p  **  are  calculated  from  the  following  version 
of  equations  (B.4)  and  (B.S), 


where  c^jt  •  l/2(Cj  ♦  etc.  Equations  (b.l)  and  (B.2)  are  then 

solved  simultaneously  for  the  new  position  of  point  3,  X3**,  and  t^**, 
by  using  the  most  recent  values  u^”  and  C3*.  Similarly,  x^**  and  t^** 
are  calculated  from  equations  (B.3)  and  (B.7).  Again,  the  most 
recently  calculated  values  of  x^**,  t3*',  u^*,  and  U3**,  are  used.  With 
x^"  established,  the  interpolation  formulas,  equations  (B.8)  to  (B.IO), 
are  used  to  calculate  u^**,  c^**,  and  p^**.  Finally,  the  quantity  c^”  is 
calculated  from  equation  (B,6)  by  using  the  values  c^'*,  p^",  and  P3**. 

The  second  approximation  of  the  properties  at  points  3  and  A  is 
now  complete.  Repeating  the  procedures  in  the  second  iteration,  with 
double  primed  quantities  instead  of  single  primed  quantities  used  as 
input,  we  can  obtain  the  third  approximation.  This  iteration  process 
may  be  continued  until  the  desired  degree  of  accuracy  is  reached.  In 
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the  present  calculation^  the  coefficients  in  equations  (B.4)  and  (B.S)^ 


were  used  to  determine  the  convergence  of  the  variables  at  point  3. 

The  iteration  was  stopped  when  each  of  the  above  four  coefficients  con¬ 
verged  to  within  .005%  of  its  value  in  the  previous  iteration.  When 
this  iteration  procedure  is  used^  the  variables  and  t^  converge 
faster  than  the  properties  u^,  c^,  and  p^.  For  roost  points  in  this 
problem^  not  more  than  four  iterations  were  required  to  achieve 
the  desired  accuracy. 

A  few  other  iteration  schemes*  have  been  tested  and  conqpared  with 
the  one  just  described.  These  schemes  give  the  desired  accuracy »  but 
the  one  just  described  seems  to  be  most  rapid. 


2)  Points  on  the  Shock  Front 


Referring  to  Figure  B.2^  all  variables  are  assumed  known  at  points 
1  and  2,  and  the* variables  at  point  3  are  to  be  determined.  Points  2 
and  3  lie  on  the  shock  front.  The  finite-difference  forms  of  the  1- 
physical  and  state  characteristics  are,  from  equations  (14)  to  (17) 


*For  instance^  instead  of  solving  for  x^'*  and  t  **  from  equations  (B.3) 
and  (B.7),  we  could  solve  equations  (B?3) ,  (B.y)  and  (B,8)  simultaneously 
for  X  t  **,  and  u^**  by  using  the  previously  determined  values  x^**,  tj**, 
and  uy'.  But  this  method  involves  the  solution  of  a  quadratic  equation 
which  has  been  found  impractical. 


y 
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(B.ll) 


Fro«  the  definition  of  shock  wave  velocity,  (^^/^^Jjhock  " 


(B.13) 


Also,  the  properties  at  point  3  must  satisfy  the  strong  shock  equations, 
(11)  to  (13);  therefore, 


^3=y?it/3 

c-=/SZZ3ru, 

P-3  =7?i  A  U  3 


(B.14) 

(B.15) 

(B.16) 


From  these  six  equations,  (B.ll)  to  (B.16),  the  six  unknowns,  x^,  t^, 
U3,  C3,  P3,  and  U3,  may  be  determined.  The  solution  is  accomplished 
by  the  following  iteration  procedure. 


First  Iteration 

Equations  (B.12),  (B.14),  and  (B.16)  are  solved  simultaneously  for 
U3#  P31  and  by  neglecting  the  term  containing  t  in  equation  (B.12). 
These  values  are  the  first  approximations  and  are  designateu  by  U3',  P3*, 
and  U3*.  The  quantity  c^*  is  calculated  from  equation  (B,15),  where 
is  inserted  for  U^.  The  position  of  point  3,  x^*  and  t^*,  is  found  by 
simultaneously  solving  equations  (B.ll)  and  (B.13).  In  these  equations, 
the  previously  calculated  values  u^',  c^',  and  are  used,  completing 
the  first  approximation  of  the  variables  at  point  3. 
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Second  Iteration 


Once  the  first  approximation  is  established,  the  following 
procedure  may  be  used  to  obtain  the  second  and  subsequent  approximations. 
Equations  (B.12},  (B.14),  and  (B.16)  are  solved  for  u^**,  p^**,  and  U^** 
by  using  some  of  the  values  obtained  in  the  first  iteration,  i.e., 

= jh 

The  quantity  c^**  is  not  calculated  from  eqi^ation  (B.15)  with  inserted 
for  U^.  The  position,  and  t^**,  is  found  by  simultaneously  solving 
equations  (B.ll)  and  (B.13).  Again,  the  most  recent  values  u  **,  c  **, 
and  are  used.  The  second  approximation  of  the  properties  at  point  5 
is  now  con^lete.  By  repeating  the  above  procedures,  with  double  primed 
quantities  replacing  single  primed  quantities,  we  can  obtain  the  third 
approximation.  This  iteration  process  may  be  continued  until  the  desired 
accuracy  is  obtained.  In  the  present  calculation,  the  iteration  was 
stopped  when  each  of  the  coefficients  in  equation  (B.12),  (C|3/p^3) 
and  ('ii3C23/Xj3) ,  converged  to  within  .005%  of  its  value  in  the  previous 
iteration. 


X 

Figure  B.3  Physical  Plane 


1 
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Referring  to  Figure  all  variables  at  points  1  and  2  are 
asstmed  knoim,  and  the  variables  at  point  3  are  to  be  detemined. 
Points  1  and  3  are  on  the  back  boundary.  Point  2  is  an  interior 
pointi  and  points  1  and  2  are  on  the  sane  I*characteristic,  Point  A 
is  located  at  the  intersection  of  the  Ill-characteristic  passing 
through  point  3  and  the  I-characteristic  joining  points  1  and  2, 

The  properties  at  point  A  will  be  determined  by  a  linear  interpolation 
of  the  properties  at  points  1  and  2,  Hie  finite-difference  forms  of 
the  physical  characteristics  are^  from  equations  (14)  and  (16) 


along  11 

(B.17) 

along  III 

(B.18) 

The  corresponding  state  characteristics  are 

along  II 

-3)  ^  (t,  tj) 

(B.19) 

along  III 

{B.20) 

Since  the  I-characteristic  is  assumed  to  be  a  straight  line  between 
points  1  and  2,  we  have 

Point  3  must  lie  on  the  back  boundary;  therefore,  it  must  satisfy  the 
equation  of  the  back  boundary  curve  presented  in  section  III. 

X3=K(^-t3^  (».«) 

where  k  ■  constant  (the  position  ratio).  Also,  as  discussed  in  section  III, 
the  particle  velocity  at  points  on  the  back  boundary  is  prescribed  by  the 
exact  solMtion;  therefore. 
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(B.23) 


where  F  »  constant  (the  velocity  ratio  •  The  linear  interpolation 

formulas f  which  are  used  to  determine  the  properties  at  point  A,  are 


(B.24) 

(B,25) 

(B.26) 


From  these  ten  equations,  (B«17)  to  (B.26),  the  ten  unknowns,  X3,  t3, 
^3*  ^3*  ^3*  '^A*  ^A*  ^A*  determined.  The  solution 
is  accomplished  by  the  following  iteration  procedure. 


First  Iteration 

Equations  (B.17)  and  (B.22)  are  solved  simultaneously  for  and  t^ 

by  using  u^  for  and  c^  for  These  values  are  then  designated  by 

X3*  and  t^*,  respectively.  The  first  approximation  of  the  particle  velocity 

u^*  is  determined  by  using  t^*  in  equation  (B.23}  and  solving  for  u^. 

The  quantity  p  '  is  determined  from  equation  (B»19)  by  neglecting  the 

term  containing  t  and  using  the  properties  at  point  2  in  place  of  the 

average  properties.  Next,  x^'  and  t^*  are  found  by  solving  equation  (B.18) 

and  (B.21)  with  u^^  replaced  by  u^*.  The  properties  u^*,  c^',  and  p^* 

are  determined  from  the  interpolation  formula.;,  equations  (B.24),  (B.25), 

and  (B.26),  and  the  calculated  position  x.*.  The  sound  speed  c  *  may  now 

^  3 

be  calculated  from  equation  (B.20),  The  first  approximation  of  the 
properties  at  points  3  and  A  is  now  coaq^lete. 
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Second  Iteration 


The  following  procedure  may  be  used  to  calculate  the  second  and 
sti>sequent  approxinations. 

Equations  (B.17)  and  (B.22)  are  solved  for  x^"  and  t^"  by  using 
the  average  properties  u^^*  and  particle  velocity  u^”  is 

calculated  froa  equation  iB«23)  by  using  t^".  Next^  equation  (B«19)  is 
solved  for  p^'*  by  using  soae  of  the  values  obtained  in  the  first 
iteration^  i-c., 


The  position  of  point  x^”  and  is  detemined  from  equation  (B.13] 
and  (B.21)  with  u^„^,  used  as  the  average  particle  velocity.  The  position 
of  point  A  nay  also  be  determined  by  the  alternate  procedure  described 
in  the  footnote  of  the  interior  point  procedure.  Next,  the  interpolation 
fonmlas,  (B.24)  to  (B.26),  arc  used  to  determine  u^\  c^**,  and  p^*’. 

Then  the  sound  speed  c^"  is  calculated  from  eouation  (B.20),  completing 
the  second  approximation  of  the  properties  at  points  3  and  A.  As 
described  in  the  general  point  iteration  procedure,  this  process  may  be 
continued  until  the  desired  accuracy  is  reached.  In  the  present  calculations, 
the  iteration  was  stopped  when  the  coefficients  in  equation  (B.19), 

(C23/P23)  and  (u23^23^^23^»  converged  to  within  .005%, 
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APPENDIX  C 


Error  Analysis  of  the  Finite«Difference  Eqitmtions 

It  is  desirable  to  determine  the  nature  of  the  truncation  error 
which  is  introduced  by  replacing  the  governing  differential  character¬ 
istic  equations  (4),  (S),  and  (6)  by  the  finite- difference  equation^ 
(14)^  (15) p  and  (16) »  Since  these  equations  are  of  the  sa.^  general 
form,  only  one  typical  equation  will  be  analyzed.  This  equation  is 
the  I-state  characteristic  equation  (6) . 

where  a  represents  arc  length  along  the  I -characteristic.  When 
this  equation  is  expressed  in  finite-difference  form,  the  coefficients 
c/p  and  cu/x  may  be  approximated  by  either  c/p  and  cu/x  or  (c/p)  and 
(cu/x) ,  i.e.,  either  by  averaging  the  variables  individually,  or  by 
averaging  the  coefficients  which  are  combinations  of  the  variables. 

It  will  be  shown  that  the  types  of  error  involved  in  both  cases  are 
the  same. 

Referring  to  Figure  C.l,  the  finite-difference  characteristic 
equation  (15)  is  to  be  written  for  point  2,  in  terms  of  values  of 
the  variables  at  points  I  and  3,  both  Aa  distance  away  from  2  along 
the  I-characteristic, 


Figure  C.l  Physical  Plane 
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With  individual  variables  averaged «  eq.  (15)  becomes 


where  t  represents  the  central  difference,  i.e.,  Au  ■  ;  and 

the  barred  quantities  represent  averages  between  points  1  and  3,  e.g., 

u»=^(ai+Uj) 

The  Taylor  series  expansions  (assuming  they  exist)  for  and  Uj 
about  point  2  are, 

ijc « i- . . . 

The  other  variables  at  points  1  and  3  can  be  expanded  in  a  similar 
manner e  By  adding  and  subtracting  the  above  series  expansions,  and 
rearranging,  we  obtain 

U,=Ua-|,(-^|^d<X*+  •  -  • 

Again,  similar  expressions  for  C2,  P2#  X2,  (dp/da)2  and  (dt/do)2  can 
be  formed  from  the  corresponding  series  expansions. 

Substituting  these  expressions  into  the  exact  differential 
equation  (C.l),  which  is  valid  at  point  2,  we  obtain 
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where  etc.,  are  independent  of  Aa  and  depend  only  upon  the 

position  of  point  2.  By  cowparing  eq.  {C.3)  with  (C.2),  we  observe 
that  the  terns  Nj(Ac)^  ♦  N^CAo)**  ♦  ...  are  neglected  in  eq.  {C.2). 
Therefore,  the  error  e  introduced  by  the  replacenent  of  eq.  (C.l)  by 
{C.2)  is  of  the  order  (Ao)^,  or 

e  *  (C.4) 

With  coefficients  averaged,  equation  (IS)  becones 

where 

= HT*  ^ 

For  this  case,  the  quantities  c/p,  cu/x,  u,  p,  and  t  at  points  3  and  1 
are  expanded  into  Taylor  series  about  point  2. 

By  combining  these  two  equations,  we  have 

-  Ml}'  ^  'ii 

If  we  substitute  the  above  equation  and  similar  expressions  for 
(dp/da) 2,  (cu/x)^  and  (dt/da)^  into  the  exact  equation  (C.l),  we  obtain 
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where  K^,  iC2t  etc.,  are  independent  of  Aa.  Coopering  eq.  (C.6)  with 
eq.  (C.5),  it  is  observed  that  the  teros  Kj(Aa)^  ♦  K^CAa)**  ♦  ...  are 
neglected  in  eq.  (C.5).  Therefore,  the  replacenent  of  (C.l)  by  (C.5) 
introduces  an  error  e'  of  order  (Aa)^  , 

We  conclude  thet  for  either  of  the  two  methods  of  averaging »  the 
errors  involved  in  all  the  finite-difference  equations  are  of  the  (Aa)^- 
type.  If  h  is  a  quantity  which  is  in  some  way  proportional  to  the  arc 
length  Ao,  then  these  errors  are  of  the  h^  -type. 
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APPENDIX  D 


Siailarity  Solution  of  the  Blast  Wave  Problem 

In  this  appendix^  Sedov's  solution  for  the  wave  generated  by 
the  instantaneous  release  of  a  finite  amount  of  energy  is  sunstarized 
(see  ref.  (13),  chap.  IV), 

From  dimensional  considerations,  the  dependent  variables,  u,  p, 
and  p,  may  be  expressed  as  follows: 


u  v(x) 


(0.1) 

(D.2) 

(D.3) 


where 


X 


(D.4) 


is  a  nondimensional  variable  and  a,  K,  s,  b^  m,  and  6  are  constants. 

When  equations  (D.l),  (D.2),  and  (D.3)  are  substituted  into  the  equations 
governing  the  one-dimensional  flow  of  a  perfect  gas  without  heat  transfer 
and  friction,  equations  (1),  (2),  and  (3),  the  following  equations  result. 


v(v-ij^v%'+|/-yv;z 


(D.6) 


(U.7) 


where 
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These  equations  are  ordinary  differential  equations  and  motions 
corresponding  to  them  are  called  self-similar «  In  these  equations^ 
a  new  variable  Z  »  y(P(X)/R(X))  has  been  introduced  and  the  expression 
for  X,  equation  (D.4),  has  been  used. 

In  the  blast  wave  problem^  a  «  and  b  «  £/p^«  where  p^  is 

the  initial  density  of  the  gas  outside  the  wave  and  E  is  a  constant 
which  is  proportional  to  the  energy  released.  The  exponents  6  and 
m  are  determined  from  the  dimensions  of 


-Vf  2 


Also,  the  constants  K  «  -3  and  S  ■  0  are  determined  from  the  dimension 

of  p..  At  the  shock  front,  x  »  x  and  X  *  X  ;  therefore,  from  equation 
'  s  s 

(D.4),  the  equation  of  the  shock  front  becomes 


(D.8) 


where  X  has  been  arbitrarily  set  equal  to  unity  at  the  shock  front, 
s 

At  the  shock  front,  the  variables  u,  p,  and  p  must  satisfy  the 
strong  shock  equations,  which  are,  in  terms  of  the  new  variables. 


These  strong  shock  boundary  conditions  and  the  stipulation  that  the 
energy  within  the  wave  zone  remains  constant  determine  a  unique  solution. 
With  the  above  information  established,  the  properties  at  the  shock 
front  become 

(D.IO) 


(Dell) 

(D.12) 
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Sedov  obtained  a  solution  of  the  differential  equations  (0.5) , 
(D.6)»  and  (D.7)  which  satisfies  the  boundary  conditions  (D.9)«  This 
solution »  when  transforaed  back  to  diaensional  variables,  aay  be 
expressed  in  teras  of  a  paraaeter  V  as  follows,  (see  ref«  13,  page  219). 


::Jy±2UJt!J^V4- 


(0.13) 

(0.14) 


(0,15) 

(0.16) 


" 0>iri ”  ®‘*]  ^2  "  srrSU' 

Ifr^hT 

The  parameter  V  has  the  range 

"r-v+liT"^  ^  TVJTzj^jj 


for  1  <  Y  <  7. 

The  complete  solution  of  the  blast  wave  problem  is  given  by 
equations  (D.IO)  through  (D,16)  together  with  the  equation  of  the 
shock  front  (D*8). 
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TABUE  11« 


x/x^ 

1. 00000000 
.99000008 
.9B000006 
.9700000S 
.96000003 

.94999993 

.94000007 

.93000011 

.92000017 

.91000011 

.90000013 

.89000010 

.88C00013 

^87000010 

.8600001*5 

.R50C0011 

.83999997 

.83000010 

.82000017 

.81000006 

.80000021 

.79000022 

.7B00001S 

.77000018 

.76000011 

.75000007 

.74000015 

.73000019 

.72000025 

.71000019 

.70000014 

.6900001? 

.60000009 

.67000008 

.66000037 

.65000005 

.64000030 

.63000038 

.62000015 

.61000030 

.60000031 


BLAST  NAVE  PROPERTIES  FOR  PLANE.  CYLINDRICAL.  AND  SPHERICAL 
NAVES  WITH  y  -  7/5  and  5/3 


a.  Plane  Nave  (v*l),  y  ■  7/5 


u/u. 

0/P3 

p/p. 

1.00000000 

1. 00000000 

1.00000000 

1.000000 

.98507586 

.92860073 

.95681441 

1.015077 

.97030586 

.86384121 

.91702569 

1.030323 

.95569410 

.80496880 

.88031874 

1.045756 

.94124416 

.75132872 

.84641294 

1.061392 

.92695898 

.70234947 

.81505749 

1.077252 

.91284169 

.65753374 

.78602963 

1.093353 

.89089401 

.61644331 

.75912832 

1.109713 

.88511779 

.5F869431 

.73417377 

1.126353 

.87151407 

.54394807 

.71 100374 

1.143292 

.85808381 

.51190647 

.68947255 

1.160548 

.84482706 

.48230480 

.66944005 

1.178142 

.83174370 

.45490900 

.65081098 

1.196093 

.81883284 

.42951019 

.63345260 

1.214423 

.80609347 

.40592327 

.61727472 

1.233153 

.79352378 

.38398229 

.60218737 

1.252305 

.78112174 

.36353930 

.58810889 

1.271900 

.76888539 

.34446246 

.57496507 

1.291962 

.75681160 

.32663233 

.56268731 

1.312514 

.74489716 

.30994204 

.55121304 

1.3335P1 

. 73313927 

.29429640 

.54048565 

l.3551f:.8 

.72153383 

.27960816 

.53045204 

1.377362 

.71007714 

.26579956 

.52106416 

1.4001 30 

.69876538 

.25280027 

.51227773 

1.423521 

.68759415 

.24054636 

.5040517*t 

1.447565 

.67655934 

.22898031 

.49634069 

1.472293 

.66565660 

.21804985 

.48913387 

1.497738 

.65488133 

.20770733 

.48237509 

1.523935 

.64422905 

.  19790966 

.47604272 

1.550920 

.63369507 

.10861748 

.4701092^ 

1.578737 

.62327496 

. 1 7979508 

.46454921 

1.60741  1 

.61296423 

,17141004 

.45933918 

1.636999 

.60275827 

.16343250 

.45445/18 

1.667543 

.59265269 

.15583539 

.44988294 

1.699091 

.58264331 

.14859405 

.44559759 

1.731692 

.57272490 

.  14168517 

.44158317 

1.765403 

.56289435 

.13500878 

.43782378 

1.800279 

.55314676 

.12078528 

.43430378 

1.836384 

.54347792 

.12275718 

.43100834 

1.873784 

.53388451 

.11698885 

.42792592 

1.912546 

.52436215 

.11146524 

.42504223 

1.952748 
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TABLE  II.  a. 

(continued) 

x/x^ 

u/u^ 

.59000036 

.51490728 

.56000030 

.50551608 

.57000026 

.49618517 

.56000061 

.48691149 

.55000062 

.47769108 

.54000071 

.46852116 

.5300004? 

.45939827 

.52000038 

.45032008 

.51000059 

.44128379 

.50000069 

.43228639 

.49000142 

.42332598 

.48000121 

.41439873 

.47000134 

.40550343 

.46000207 

. 39663310 

.45000199 

.38779939 

.44000171 

.37^98584 

.43000155 

.37019584 

.42000147 

.36142753 

.41000221 

.35267990 

.40000337 

.34395096 

.39000176 

.33523643 

.38000264 

.32653950 

.37000275 

.31785596 

.36000425 

.30918647 

.35000328 

.30052648 

. 34000270 

.29.8774? 

.33000476 

.28324017 

.32000414 

.27460923 

.31000672 

.26598870 

.30000857 

.25737460 

.29000825 

.24876474 

.28001068 

.24016272 

.27000359 

.23155736 

.26001251  ■ 

.22297004 

.25001419 

.21438023 

.24001929 

.20579662 

.23001775 

.19721014 

.22002668 

.18863509 

.2100W8 

.  18004683 

.20001246 

.17146339 

p/p. 

9 

VI 

. 10617291 

.42234623 

.10109933 

.4198268? 

.09623310 

.41747380 

.09156385 

.41527757 

.08708148 

.41322885 

.08277720 

.41131921 

.07864253 

.40954047 

.07467007 

.40783517 

.07085273 

.40634611 

.06718380 

.4049164? 

.06365748 

.40358984 

.06026761 

.40236012 

.05700940 

.40122161 

.05387808 

.40015893 

.05036870 

. 39919632 

.04797728 

.  398  *0041 

.04520000 

.39747510 

.04253319 

.39671643 

.0*997365 

.39602033 

.03751804 

.39538270 

.0*51626? 

.39479963 

.03290585 

.39426782 

.03074428 

.39378367 

.02867588 

.3^334396 

.02669751 

.39294554 

.02480748 

. 39253553 

.02300405 

. 39226123 

.02128420 

.39196980 

.01964698 

.39170892 

.01808980 

.39147606 

.01661065 

.39126895 

.01520848 

.  39108557 

.01387990 

.39092372 

.01262654 

.39070182 

.01144371 

.39065776 

•01033109 

.  39054996 

.0092P595 

.39045677 

.00830043 

.  *9037681 

.00739421 

.39030857 

.00654376 

.39025078 

St 


1.994469 

2.037797 

2*O02ft2? 

2.  129644 
2.178373 

2.229123 
2. 28202! 
2.33719P 
2.394802 
2.454993 

2.517939 

2.583839 

2.652887 

2.725308 

2.8013S5 

2.881293 

2.965416 

3.054049 

3.  147544 
3.246300 

3.350793 

3.461458 

3.578875 

3.703631 

3.836461 

3.978100 
4. 129386 
4.291385 
4.465126 
4.651956 

4.053379 

5.07099? 

5.307044 

5.563203 

5.84271A 

6. 14844? 
6.484449 
6.854603 
7.26537? 
7.722500 


TABLE  II. 

BLAST  WAVE  PROPERTIES  FOR  PLANE, 

CYLINDRICAL, 

AND  SPHERICAL  ^ 

WAVES  WITH  y  - 

7/3  and  5/3 

T 

b.  Cylindrical  Wave  (vb2),  y  ■  7/S 

1 

x/x^ 

U/U^ 

P/Pj 

p/p. 

I 

U 00000000 

1. 00000000 

1.00000000 

1 .00000000 

1. 000000 

.99000023 

.98187036 

.87082494 

.92731237 

1.0J1923  _ 

.96000020 

.96415348 

.76311560 

.86460582 

l.064«?2  1 

.97000012 

.94685728 

.67256243 

.81022513 

1.097580  * 

.96000015 

.92998646 

.5958488? 

.76283978 

1.131484 

T 

.95000013 

.91354205 

.53039435 

.721BS991 

1.166217  1 

.94000010 

.89752275 

.47417499 

.68493403 

1.201863 

.9300002? 

•86192341 

.42558459 

.65280404 

1.238507 

.92000016 

.86673740 

.38334339 

.62437809 

1.276232 

.91000027 

.85195567 

.34642050 

.59915324 

1.315125 

.90000027 

.83756658 

.31397892 

.57670605 

1.355273 

.6900003? 

.82355753 

.28533698 

.55668014 

1.396766 

•68000021 

.80991407 

.25993386 

.53877216 

1.439697 

.87000035 

.79662161 

.23730774 

.52272448 

1.484159 

•66000028 

.78366368 

.21707312 

.50831513 

1.530253 

.85000040 

.77102453 

.19890982 

.49535434 

1.578083 

.64000030 

.75668710 

.18254765 

.48367726 

1.627757 

.83000041 

.74663534 

.16776004 

.47314185 

1.679389 

.62000034 

.73485240 

.15435402 

.4636238! 

1.733101 

.81000022 

.72332229 

.14216592 

.45501505 

1.7C9019  j 

.80000037 

.71202963 

.13105615 

.44722103 

1.84727P 

.79000037 

.70095867 

.12090413 

.44015816 

1.908023  t 

.78000036 

.69009484 

.11160654 

.43375317 

1.971407  1 

.77000041 

.67942411 

.  10307392 

.42794114 

2.037593 

.76000037 

.66893280 

.09522847 

.4226645? 

2.106757 

.75000063 

.65860849 

.08800280 

.41787236 

2.179081 

.74000034 

.64843811 

.08133709 

.41351864 

2.2*^4775 

.73000047 

.63841113 

.07518002 

.40956300 

2.334044  [ 

.72000019 

.62851581 

.06948532 

.40596864 

2.417128  i 

.71000026 

.61874270 

.06421298 

.40270295 

2.504269 

.70000058 

.60908208 

.05932704 

.39973646 

2.595736  | 

.69000115 

.59952492 

a05479554 

.39704233 

2.691818 

.68000101 

.59006189 

.05058956 

.39459632 

2.792839 

.67000119  ' 

.5806861? 

.04668403 

.39237678 

2.879128  1 

.66000132 

.5713899? 

.04305583 

.39036383 

3.011055  • 

.65000129 

.56216638 

.03968425 

.38853941 

3.129020  1 

.64000164 

.55300965 

.03655071 

.38688719 

3.253450  1 

.63000149 

.54391306 

.03363795 

.38539203 

3.384827 

.62000192 

.53487216 

.03093081 

.38404040 

3.523649 

.61000184 

.52588090 

.02841485 

.38281962 

3.670492  1 

.60000282 

.51693607 

.02607750 

.38171045 

3.825943  * 

52 

TABLE  II.  b. 


x/x^ 

.*59000272 

.S8000315 

.57000373 

.56000345 

.55000444 

.54000425 

.53000744 

.52000608 

.51000749 

.50000719 

.49000826 
.4800101  I 
.47000600 
.46000808 
.45001258 

.44001434 

.43001878 

.42002251 

.41001927 

.40001952 

.39002127 

.38004135 

.37003513 

.36003655 

.35005084 

.34007201 

.33009187 

.32012595 

.31018088 

.30013331 

.29020790 

.28C25050 

.2702619^ 

.26046833 

.25049917 

.24053268 

.23059844 

.22101686 

•21141195 

.20210235 


(continued) 

P/Pj 

.50803102 

.02390626 

•49916457 

.02109082 

.49033143 

.02002066 

.48152783 

.01828630 

.47275259 

.01667937 

.46400074 

.01519136 

.45527377 

.01381533 

•44656241 

.01254309 

.43787094 

.01136897 

.42919351 

.01028617 

.42053110 

.00928916 

.41188158 

.00837233 

.40323773 

.00752990 

.39460875 

.00675797 

.3H599023 

.00605150 

. 37737669 

.00540576 

.36877186 

.00401692 

.36017201 

.00428085 

.35157103 

.00379355 

.34297723 

.0033520? 

.33438036 

.00295281 

.32581833 

.00259333 

•31722039 

.00226903 

. 30864727 

.00197830 

.30007912 

.00171853 

.29151849 

D014P699 

.28295811 

.00128112 

.27441106 

.00109897 

.26588283 

.00093849 

.25726750 

.00079598 

.24875756 

.00067276 

.24022072 

.00056497 

.23165757 

.00047120 

.22326191 

.00039179 

.21471602 

.00032233 

.20617263 

.00026310 

.19765707 

.00021307 

.18944392 

.00017233 

.18121087 

.00013800 

.17323103 

.00011018 

p/p,  c/c. 


.38072008 

3.990712 

.37983320 

4.165483 

.37903030 

4.351089 

.37831073 

4.548428 

.37766569 

4.758431 

.37708865 

4.982224 

.37657364 

5.??0887 

.37611445 

5.475927 

. 37570614 

5.748618 

.37534364 

6.040704 

.37502261 

6. 353899 

.37473898 

6.690230 

.37440878 

7.052201 

.37426890 

7.441898 

.37407615 

7.862273 

.37390751 

8.316748 

. 37376046 

8.808695 

. 37^63264 

9.342368 

.37352181 

9.922814 

.37342614 

10.554770 

.37334385 

11.244401 

.37327339 

11.997327 

. 37321309 

12.625022 

.37316188 

13.734168 

.37311059 

14.734785 

.37308214 

15.839694 

.37305156 

17.064317 

.37302609 

18.423639 

.37300503 

19.936156 

.37298745 

21.646845 

. 37297323 

23.545504 

.37296161 

25.693058 

.37295222 

28.133270 

.37294482 

30.052869 

. 37293882 

34.014766 

.37293407 

37.648559 

.37293043 

41.835300 

.37292767 

46.517937 

.37292555 

51.982875 

.  37  292397 

58.177597 

S3 


TABLE  ll.  BLAST  WAVE  PROPERTIES  FOR  PLANE,  CYLINDRICAL,  AND 
WAVES  WITH  Y  7/5  and  5/3 


c«  Spharical  Wava  (v«3),  y  *  7/5 


x/x^ 

u/u^ 

p/pj 

P/Pj 

1.00000000 

l.OOOCOOOO 

1.00000000 

1.00000000 

.99000020 

.97872610 

.81781524 

.89971703 

.98000016 

.95824629 

.67765852 

.81849516 

.97000023 

.93856766 

.56797341 

.75197047 

.96000031 

.91968545 

.48080899 

.69694604 

.95000024 

.90158501 

.41058013 

.65103994 

.94000019 

.88424384 

.35328933 

.61244954 

.93000016 

.86763231 

.30602359 

.57978968 

.92000028 

.85171574 

.26662837 

.55198323 

.91000032 

.83645521 

.23348568 

.52818167 

.90000031 

.82181004 

.20536700 

.50771146 

.89000036 

.80773849 

.18132721 

.49003170 

.88000052 

.79419876 

.16063148 

.47470484 

.87000048 

.78114922 

.14270158 

.46137315 

.86000043 

.76855028 

.12707988 

.44974330 

.85000054 

.75636395 

.11339969 

.43957251 

.84000052 

.74455379 

.10136447 

.43065770 

.83000062 

.73308609 

.09073326 

.42282925 

.82000066 

.72192898 

.08130799 

.41594376 

.81000076 

.71105318 

.07292514 

.40987984 

.80000076 

.70043129 

.06544827 

.40453367 

.79000072 

.69003843 

.05876321 

.39901652 

.78000095 

.67985193 

.05277360 

.39565206 

.77000115 

.66985054 

.04739735 

.39197405 

.76000103 

.66001482 

.04256423 

.38872502 

.75000132 

.65032811 

.03821437 

.38585510 

.74000120 

.64077385 

.03429527 

.38332042 

.73000157 

.63133853 

.03076198 

.38108273 

.72000170 

.62200855 

.02757462 

.37910813 

.71000153 

.61277216 

.02469836 

.37736680 

.70000246 

.60362011 

.02210287 

.37583259 

.69000334 

.59454187 

.01976056 

.37448200 

.68000320 

.50552805 

.01764701 

.37329420 

.67000259  • 

.57657152 

.01574068 

.37225086 

.66000479 

.56766836 

.01402270 

.37133587 

.6500031 9 

.55880658 

.01247447 

.37053408 

.64000410 

.54998634 

.01108120 

.36983305 

.63000660 

.54120193 

.00982838 

.36922099 

.62000647 

.53244535 

.00870257 

.36868736 

.61000535 

.52371418 

.00769224 

.36822301 

.60000740 

.51500870 

.00678706 

.36781995 
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SPHERICAL 


c/c^ 

1. 000000 
1.04«a79 
1 .099012 
1.150631 
1.203963 

1.259229 

1.316648 

1.376441 

1.438830 

1.504047 

1.572328 

1.643919 

1.719081 

1.798091 

1.881239 

1.968835 

2.061214 

2.158734 

2.261781 

2.370771 

2.406155 

2.608420 

2.738093 

2.875753 

3.022030 

3.  177596 
3.343211 
3.519674 
3.707889 
3*908838 

4. 123564 
4.353271 
4.599282 
4.863019 
5.145973 

5.450081 

5.777094 

6.129174 

6.508863 

6.918770 

7.361676 


TABLE  II.  c*  (continued) 


.59000B39 

•580014P0 

.57001584 

.56001335 

.55001672 

.54002796 

.53003275 

.52002840 

.51002923 

.50004599 

.49003355 

.48009163 

.47006348 

.46010146 

.45011067 

.44018667 

.43015107 

.42019308 

.41030184 

.40030712 

.39047057 
.38050044 
.3707871  1 
.360M430 
.35088826 

.34116663 

.33146455 

.32133957 

.31262625 

.30279608 

.29354415 

.28450148 

.27626264 

.26655389 

.25903206 

.25130555 

.24080536 

.23U9252 

.22073661 

.20190170 


U/Uj 

.50632217 

.49765766 

.48900326 

.48035874 

.47173035 

.46311823 

.45450867 

.44589817 

.43729803 

.42871658 

.42011425 

.41157605 

.40296679 

.39441679 

.38584417 

.37733058 

.36872267 

.36018251 

.35170054 

.34313061 

.33469694 

.32614924 

.31762215 

.30920425 

.30076393 

.29243045 

.28411388 

.27543493 

.26796611 

.25954004 

.25160967 

.24385870 

.23679676 

.22847491 

.22202759 

.21540483 
.20640464 
.19842220 
.1892  0283 
.17305860 


p/p^ 

.00597684 

.00525316 

.00460709 

.00403143 

.00351990 

.00306626 

.00266420 

.00230854 

.00199508 

.00171971 

.00147732 

.00126657 

.00108094 

.00092039 

.00078058 

.00066033 

.00055540 

.00046589 

.00038965 

.00032384 

.00026870 

.00022132 

.00018229 

.0001^4833 

.00012053 

.00009762 

.00007863 

.00006230 

.00005070 

.00003989 

.00003161 

.00002500 

.00002005 

.00001533 

.00001237 

.00000985 

.00000715 

.00000532 

.00000372 

.00000190 


p/p, 

.36747065 
.36716078 
.36690825 
. 3666H391 
.36649147 

.36632683 

.36618617 

.36606627 

.36596456 

.36587863 

.36580593 

.36574522 

.36569391 

.36565139 

.36561595 

.36550674 
. 3655624] 
.36554262 
. 36552643 
.36551325 

.36550271 

.36549408 

.36548730 

.36548174 

.36547743 

.36547408 

.36547145 

.36546932 

.36546787 

.36546664 

.36546573 
.36546504 
.36546458 
.36546414 
. 36546390 

.36546369 

.36546348 

.36546337 

.36546327 

.36546317 


C/C^ 

7.841068 

8.360305 

8.924111 

9.537092 

10.203909 

10.930227 

11.723759 

12.592464 

13.543761 

14.586153 

15.735750 

16.993151 

18.393160 

19.931790 

21.642277 

23.529600 

25.655243 

28.010654 

30.628093 

33.595738 

36.881095 

40.637672 

44.775833 

49.637383 

55.065425 

61.184002 

60.174751 

76.585588 

84.901924 

95.708296 

107.519620 

120.905540 

134.991110 

154.372310 

171.865060 

192.533140 

225.951510 

261.968880 

313.137030 

437.510460 
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TABLE  !I.  BLAST  WAVE  PROPERTIES  FOR  PLANE,  aLINDRICAL,  AND  SWIERICAL 
WAVES  WITH  Y  «  7/5  «nd  5/3 


d.  Plane  Wave  (v>l),  y  •  5/3 


x/x. 

u/u. 

P/Pg 

p/p. 

c/c. 

1« 00000000 

1. 00000000 

1. 00000000 

1.00000000 

1. 000000 

^99000004 

.98505072 

.95624785 

.96593220 

1.005051 

«9B000000 

.97020531 

.91488338 

.93365708 

1.010208 

•97000004 

.95546789 

.87575597 

.90307468 

1.015477 

•96000010 

.94084226 

.83872493 

.87409084 

1.020865 

•95000006 

.92633199 

.00365084 

.84661692 

1.026378 

•94000001 

.91194081 

.77043567 

.82057006 

1.032023 

•93000004 

.89767237 

.73894217 

.79537791 

1.037807 

.9?000005 

.88352987 

.70907195 

.77245214 

1.043735 

•91000002 

.86951638 

.68072551 

.75023869 

1.049817 

.90000001 

.85563498 

.65381094 

.72916842 

1.056058 

.89000005 

•84188842 

.62824164 

.70918026 

1.062465 

.88000000 

.82827901 

.60393647 

.69021671 

1.069047 

;8700006? 

.81480920 

.58082015 

.67222393 

1.075811 

•86000001 

.80148084 

.55882175 

.65515035 

1.082764 

.85000004 

.78829569 

.53787520 

.63894953 

1.089914 

.84000000 

.77525502 

.51791799 

.62357420 

1.097270 

.83000001 

.76236011 

.49889235 

.60898232 

1.104538 

.82000002 

.74961171 

.48074356 

.59513321 

1.112628 

.80999997 

.73701023 

.46342020 

.58198839 

1.120649 

.79999998 

.72455605 

.44687477 

.56951180 

1.128907 

.78999993 

.71224891 

.43106192 

.55766934 

1.137414 

.77999997 

.70008862 

.41593945 

.54642029 

1.146176 

.76999991 

.68807420 

.40146760 

.53*^75806 

1.  155205 

.75999991 

.67620513 

.38760931 

.52562971 

1.164509 

.74999990 

.66447982 

.37432944 

.51601572 

1.174098 

.73999988 

.65289693 

.36159506 

.50685992 

1.183983 

.73000001 

.64145490 

.34937559 

.49822793 

1.194174 

.71999990 

.63015135 

.33764170 

.49000622 

1.20468? 

.70999996 

.61898438 

.32636636 

.48220280 

1.215520 

.69999984 

.60795186 

.31552303 

.47479666 

1.226698 

.68999990 

.59705104 

.30509038 

.46776821 

1.230229 

.67999995 

.58627923 

.29504327 

.46109864 

1.250126 

.66999999 

.57563358 

.20536123 

.45477007 

1.262404 

.65999985 

.56511094 

.27602415 

.44876557 

1.275076 

.64999998 

.55470874 

.26701384 

.44306958 

1.283157 

.63999989 

.54442325 

.25831223 

.43766672 

1.301665 

.62999985 

.53425153 

.24990297 

.43254274 

1.315615 

.61999985 

.52419024 

.24177067 

.42768425 

1.330025 

.60999984 

.51423596 

.23390066 

.42307830 

1.344914 

.59999976 

.50430522 

.22627920 

.41071273 

1.360303 
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TABLE  11*  d. 

(continued) 

x/x^ 

u/u^ 

p/pj 

P/Pg 

.5e999‘^^'> 

.49463495 

•21889369 

.41457608 

1.376212 

•SCOOOUOl 

.48498122 

•21173179 

.41065739 

1.392665 

.57000001 

.47542076 

.20478224 

.40694625 

1.409686 

.56000005 

.46595017 

.19803451 

.40343289 

1.427299 

.54999997 

.45656591 

.19147848 

.40010791 

1.445534 

.54000000 

.44726477 

.18510496 

.39696258 

1.464420 

.53000005 

.43804331 

.17890517 

.39398845 

1.483987 

.52000005 

.42889816 

.17287075 

. 39117741 

1.504271 

.51000004 

.41982607 

.16699411 

.^8852203 

1 .52  5306 

.49999999 

.41082478 

.16126797 

.38601499 

5.547134 

.49000026 

.40188847 

.15568580 

.38364957 

1.569794 

.48000027 

.39301649 

.15024095 

.38141905 

1.593334 

.47000000 

.38420481 

. 14492750 

.37911718 

1.617803 

.46000021 

.37545111 

.13974033 

.37n>38l7 

1.643253 

.45000020 

.366751^0 

•13467396 

.37547631 

1.669742 

.44000026 

.35810458 

.12972373 

.37372617 

1.697333 

.43000014 

.34950^  2  3 

.12498303 

.37206250 

1.726094 

.42000019 

.34095443 

. 12015390 

.37054051 

1.75609P 

.41000013 

.33244640 

.11552638 

.36909537 

1.7B7428 

.40000044 

.32397999 

.11099911 

.36774262 

1.820171 

.35999994 

.3155518? 

.10656829 

.36647776 

1.854427 

.38000004 

.30716069 

.10223150 

.36529679 

1.690299 

. 47000047 

.29880413 

.07798587 

. 36419571 

1.927905 

.36000020 

.29047908 

.09382843 

. 36317053 

1.967379 

.35000024 

.28218428 

.08975730 

.36221767 

2.009861 

.34000022 

.27391 743 

.08577027 

•361 33355 

2.052512 

.32999971 

.26567621 

.08186531 

.36051473 

2.098512 

.31999965 

.25745954 

.078041  13 

.35975787 

2.147056 

.31000081 

.24926627 

.07429656 

.35905993 

2.198361 

.30000144 

.24109326 

.07062970 

.35841767 

2.252687 

.29000062 

.23293815 

.06703917 

.35782806 

2.310323 

.280001C9 

.22480163 

.06352510 

.35728839 

2.371574 

.27000146 

.21668109 

.06008631 

.35679579 

2.436812 

.26000204 

.20857539 

.05672241 

.35634753 

2.S0645? 

.25000177 

.20048237 

.05343276 

.35594105 

2.580983 

.24000181 

.19240171 

.05021761 

.35557377 

2.660950 

.23000253 

.18433254 

.04707715 

.35524333 

2.746995 

.22000282 

.17627286 

.04401127 

.35494726 

2.939879 

.21000324 

.16822210 

.04102060 

.35468327 

2.940487 

.20000331 

.  16017889 

.03810565 

.35444912 

3.049876 

I 
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TABLE  II.  BLAST  WAVE  PROPERTIES  FOR  PLANE,  CYLINDRICAL,  AND  SPHERICAL 
HAVES  WITH  Y  -  7/5  and  5/3 

e.  Cylindrical  Nave  (v>2),  y  ■  5/3 


xfx^ 

u/u^ 

p/p. 

p/p. 

c/c. 

1. 00000000 

1.00000000 

1. 00000000 

1. 00000000 

1.000000 

.9900000? 

.98260494 

.92166523 

.94511996 

1.012644 

.98000006 

.96542669 

.85100432 

.89511844 

1.025591 

.97000009 

.94047524 

.78713225 

.04950298 

1.038863 

.96000007 

.93175984 

.72927755 

.80783703 

1.052483 

.9500000? 

.91528904 

.67676787 

.76973429 

1.066474 

.94000000 

.89907062 

.62901549 

.73485168 

1.080859 

.93000006 

.88311131 

.58550493 

.70288336 

1.095661 

.92000002 

.86741644 

.54578281 

.67355607 

1.110905 

.91000004 

.85199084 

.50945142 

.64662639 

1.126614 

.90000004 

.83683774 

.47615900 

.62187548 

1.142814 

.89000003 

.82195941 

.44559531 

.59910741 

1.159530 

.88000006 

.80735703 

.41748588 

.57814596 

1.176787 

.86999999 

.79303034 

.39158689 

.55383226 

1.194611 

.85999995 

.77897846 

.36768269 

.54102353 

1.213029 

.84999999 

.76519925 

.34558121 

.52459068 

1.232069 

.84000000 

.75168948 

.32511128 

.50941689 

1.251758 

.83000007 

.73844518 

.30612033 

.49539640 

1.272126 

.81999999 

.72546130 

.28847220 

.48243349 

1.293203 

.80999998 

.71273252 

.27204546 

.47044149 

1.31501Q 

.80000000 

.70025258 

.25673112 

.45934143 

1.337607 

.78999995 

.68801460 

.24243153 

.44906157 

1.361000 

.73000017 

.67601175 

.22905966 

.43953704 

1.385234 

.77000009 

.66423576 

.71653616 

.43070796 

1.410347 

.76000019 

.65267930 

.20479105 

.42252064 

1.436377 

.75000005 

.64133377 

.19376020 

.41492532 

1.463365 

.74000010 

.63019136 

.  18338679 

.40787716 

1.491353 

.73000008 

.61924347 

.17361896 

.40133493 

1.520389 

.72000010 

.60848179 

.16441005 

.39526082 

1.550521 

.71000016 

.59789798 

.15571789 

.38962028 

1.581799 

.70000033 

.58748382 

.  14750441 

.38438176 

1.614200 

.69000029 

.57723078 

.13973474 

•37951601 

1.648021 

.69000021 

.56713100 

.1 3237764 

. 37499632 

1.6H3085 

.67000017  ' 

.55717667 

.12540470 

.37079827 

1.719538 

.66000038 

.54736033 

.11879019 

.36689934 

1.757450 

.65000017 

.53767383 

.11251007 

.3632/830 

1.796900 

.64000014 

.5281105? 

.10654330 

.35991631 

1.837966 

.63000018 

.51866326 

.10007018 

.35679537 

1.880737 

.62000043 

.50932541 

.09547291 

.35389918 

1.925305 

.61000062 

.5000902? 

.03033503 

.35121246 

1.971772 

.60000025 

.49095098 

.08544131 

.34872098 

2.020249 
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TABLE  II,  e,  (continued) 

x/Xs  u/u^ 

.59000029  .48190255 
.58000017  .4729^869 
.57000016  .46405415 
.56000046  .45524386 
.55000088  .44650267 

.54000036  .43782492 
.53000090  .42920783 
.52000070  .42064560 
.51000123  .41213543 
.50000133  .40367251 

.49000112  .39525333 
.48000021  .35687418 
.47000121  .37853403 
.46000109  .37022733 
.45000191  .36195297 

.44000148  .-537:  48 
.43000259  .  34«;.4r?68 
.42000121  .3372  ^096 
.4100042?  .329H977 
.40000250  .32096464 

.39000401  .31?«30?1 
. 38000526  .3047UB7 
.37000562  .29660753 
. 36000438  .28851511 
.35000671  .28043744 

.34000493  .27236/03 
.33000831  .26431021 
.32001099  .25626119 
.31000811  .24821510 
.30000598  .24017854 

.29001115  .23214874 
.28000688  .22411877 
.27001044  .21609942 
.26001364  .20808353 
.25001894  .20007255 

.24002711  .19206663 
.23003225  .18406061 
.22001616  .17603958 
.21001715  .16803388 
.20007142  .16007221 


p/p. 

p/p. 

c/c^ 

.08077850 

.34641181 

2.070849 

.07633376 

.34427269 

2.123698 

.07209565 

.34229230 

2.178934 

.06805356 

.34046031 

2.236700 

.06419762 

.33076663 

2.297157 

.06051838 

.33720202 

2.360485 

.05700820 

.33575018 

2.426859 

.05365870 

.33442688 

2.496493 

.05046317 

.33320071 

2.569600 

.04741456 

. 33207261 

2.646431 

.04450670 

.33103586 

2.727249 

.04173370 

. 33008439 

2.812347 

.03909083 

.32921251 

2.902021 

.03657221 

.32841457 

2.996646 

.04417356 

.32768542 

3.096587 

.03188990 

.32702028 

3.202290 

.02971761 

.32641468 

3.314193 

.02765183 

.  .42586406 

3.432861 

.O;  >69014 

.32536488 

3.558787 

.02382713 

.32491281 

3.692728 

.02206090 

.32450461 

3.835294 

.02038744 

.32413682 

3.987335 

.01880351 

.32380626 

4.149761 

.01730599 

.32350992 

4.  32  4590 

.01580270 

.32324520 

4.509903 

.01455969 

.32300917 

4.710112 

.01330545 

. 32279965 

4.925515 

.01212651 

.32261421 

5.157910 

.01101974 

.32245057 

5.409359 

.00998364 

.32230691 

5.681856 

.00901540 

.32218128 

5.978021 

.00811186 

.32207181 

6.301093 

.00727188 

.32197698 

6.654094 

.00649230 

.32189529 

7.C41379 

.00577096 

.32102526 

7.46768? 

.00510553 

.32176562 

7.93869? 

.00449333 

.32171509 

8.461580 

.00393112 

.32167258 

9.045333 

.00341882 

.32163717 

9.699401 

.00295551 

. 42160810 

10.431497 
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TABLE  II.  BLAST  WAVE  PROPERTIES  FORPLANE,  CYLINDRICAL,  AND  SPHERICAL 
NAVES  WI1H  y  •  ?/5  and  5/3 

f,  Spharical  Wave  (v«3),  y  ■  5/3 


u/u^ 

p/pj 

p/p. 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

l.OOOOOO 

•99000010 

.98017889 

.88868587 

.92510274 

1.020283 

•9800000S 

.96072915 

.79278986 

.85941096 

1.041169 

.97000007 

.94167081 

.70980266 

.80161795 

1.062710 

.96000004 

.92302074 

.63767169 

.75062777 

1.084960 

.95000008 

.90479347 

.57471311 

.70551948 

1.107972 

.94000008 

.88700008 

.51953623 

.66551378 

1.131802 

.93000007 

.86964864 

.47098843 

.62994924 

1.156505 

•92000010 

.85274409 

.42811057 

.59026253 

1.182137 

.91000003 

•83628788 

•39010008 

.5699/1 14 

1.208755 

.90000005 

.82027884 

.35628416 

.54466125 

1.236417 

.89000013 

.80471265 

.32609564 

.52197614 

1.265181 

.8800C jll 

.78958210 

.29905463 

.50160716 

1.295109 

^87000009 

.77487813 

.27475485 

.48328734 

1.326264 

.86000006 

.76058942 

.25285034 

.46678460 

1.358709 

.85000020 

.74670322 

•23304612 

.45189683 

1.392511 

.84000015 

.73320466 

.21508865 

.43844678 

1.427741 

.8300001? 

.72007862 

.19876139 

.42628020 

1.464473 

.82000013 

.70730884 

.18307740 

.41526138 

1.502783 

.81000018 

.69487865 

.17027558 

.40527100 

1.542753 

.80000020 

.68277089 

.15781587 

.39620385 

1.584470 

.79000024 

.67096855 

.14637710 

.38796702 

1.620024 

.78000031 

.65945459 

.13585372 

.3804783? 

1.673513 

.77000025 

.64821210 

.12615344 

.37366485 

1.721042 

.76000020 

.63722480 

.11719567 

.36746182 

1.770721 

.75000020 

.62647690 

.10890978 

•36181165 

1.822669 

.74000025 

.61595267 

.10123320 

.35666263 

1.077013 

.73000041 

.60563763 

.09411123 

.35196093 

1.933880 

.72000020 

.59551704 

.08749490 

.34760897 

1.993444 

.71000058 

.58557006 

.08134153 

.34378601 

2.055833 

.70000028 

.57580676 

.07561210 

.34022620 

2.121232 

•69000066 

.56619216 

.07027315 

.33698003 

2.189814 

.68000064 

•55672163 

.06529344 

.33401982 

2.261785 

•67000068 

.54738468 

.06064579 

.33132122 

2.337352 

.66000074 

.53817101 

•05630546 

.32886178 

2.416748 

.65000052 

.52907066 

.05225007 

.32662122 

2.500222 

•64000082 

.52007530 

.04845981 

.32456124 

2.50r.O39 

•63000073 

.51117564 

.04491603 

.32272474 

2.680497 

.62000100 

.50236441 

.04160239 

.32103651 

2.777909 

.61000092 

.49363361 

.03850352 

.31950241 

2.880625 

.60000139 

.48497713 

.03560581 

•31810961 

2.909014 
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(continued) 


TABLE  II.  f. 


x/x^ 

•59000108 

.50000109 

•57000139 

•56000187 

•55000707 

.54000755 

.53000300 

.52()007«;'6 

.51000361) 

.50000413 

.49000558 

.48000500 

•47000635 

i46000640 

.45000758 

.44001103 

.43000790 

.42000973 

.41001137 

.40001846 

.39001779 
••  38002020 
.37002481 
.36003052 
.35002743 

.34002149 
.33003543 
.32006 i %8 
.31004889 
.30004625 

.29006541 

.28008761 

.27010540 

.26018978 

.25021880 

.24024027 

.23025986 

.22038949 

.21058596 

.20062742 


u/u^ 

p/pj 

p/p. 

c/c^ 

.47638741 

.03289616 

.31684624 

3.103499 

.46785946 

.03036319 

.31570148 

3.224516 

.45938773 

.02799607 

.31466540 

3.352555 

.45096704 

.0257P479 

.31372883 

3.488151 

.44259225 

.02371996 

.31286321 

3.631901 

.43425951 

.02179310 

.31212092 

3.78444? 

.42596453 

.01999610 

.31143470 

3.946488 

.41770323 

.01832136 

.31081793 

4.110832 

.40947321 

•01676204 

.31026461 

4.302318 

.40127068 

.01531139 

.30976904 

4.497919 

.39309365 

.01396332 

.30932608 

4.706671 

.38493709 

.01271156 

.30893084 

4.929820 

.37680187 

.01155107 

.30857904 

5.168589 

.36868322 

.01047620 

.30826653 

5.424518 

.36058111 

.00948223 

.30798963 

5.699184 

.35249470 

.00856450 

.30774498 

5.994380 

.34441521 

.00771775 

.30752914 

6.312446 

.33635053 

•00693668 

.30733945 

6.055345 

.32829518 

.00622277 

.30717321 

7.025864 

.32025252 

.00556648 

.30702800 

7.426746 

.31221086 

.00496526 

.30690141 

7.861914 

.30417799 

.00441626 

.30679161 

8.334781 

.29615235 

.00391601 

•30669660 

8.849772 

.28813232 

.00346123 

.30661478 

9.411986 

.28010929 

.00304842 

.30654449 

10.027885 

.27208745 

.00267496 

•30648444 

10.703968 

.26408454 

.00233883 

.30643347 

11.446393 

.25609390 

.00203691 

.30639036 

12.264529 

.24807429 

.00176527 

.306  35393 

13.173630 

.24006451 

.00152294 

.30632348 

14.182324 

.23207367 

.00130779 

.30629813 

15.303927 

.22408649 

.00111712 

.30627717 

16.557921 

.2160968? 

.00094876 

.30625995 

17.966596 

.20816126 

•00080176 

.30624599 

19.543910 

.20018208 

.00067246 

.30623459 

21.339906 

.19219740 

.00055991 

.30622543 

23.386206 

.18421167 

.00046258 

.30621015 

25.728897 

.17631430 

.00037980 

.30621245 

28.394373 

.16847070 

.00030947 

•30620806 

31.455399 

.16050328 

.00024885 

.30620456 

35.077849 
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